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Abstract 

The density matrix renormalization group (DMRG) is a numerical method that 
optimizes a variational state expressed by a tensor product. We show that the ground 
state is not fully optimized as far as we use the standard finite system algorithm, that 
uses the block structure B • »B. This is because the tensors are not improved directly. 
We overcome this problem by using the simpler block structure B • B for the final 
several sweeps in the finite iteration process. It is possible to increase the numerical 
precision of the finite system algorithm without increasing the computational effort. 



Establishment of the density matrix renormalization group (DMRG) by White |l| is one 
of the major progresses in computational condensed matter physics. DMRG enables us to 
calculate ground states of relatively large scale one-dimensional (ID) quantum systems. |^|, 
||, [|, [|, |). Two-dimensional (2D) classical systems, Jj], ||, |], |To| and ID quantum system at 
finite temperature jy], |l2], O, |l4| have also been investigated. 

Ostlund and Rommer |lq ] examined the thermodynamic limit (iV — > oo) of the infinite 
system algorithm, and they pointed out that the block state B corresponds to a product 
of position independent tensor. It should be noted that their result does not show that 
the infinite system algorithm creates a translationally invariant — position independent — 
variational state for the whole system B»»B, where "•" denotes a bare spin variable between 
the left and the right blocks. Actually, the variational state has a slight position dependence. 
For example, the bond energy (••) at the center of the antiferromagnetic S — 1/2 Heisenberg 
chain, which is calculated by the infinite system algorithm, is lower than the exact ground 
state energy per site. jl6| Such a position dependence in the variational state spoils the 
numerical efficiency of the infinite system algorithm. []l7| As we show in the following, the 
finite system algorithm does not fully improve the variational state in the same reason. The 
purpose of this letter is to remove the source of such a numerical error, and to increase the 
numerical precision in DMRG. 

Let us briefly review the construction of the variational state, which is used in the 
standard finite system algorithm. We consider the IRF model Jl8| as a reference system, 
whose transfer matrix is written as the product of local Boltzmann weights 

rp S l--- S N Jjr S l S 2 TT/ S 2 S 3 T/f/ SiY_lSjV ( 1 "\ 
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where W represents the IRF weight. The finite system algorithm, that uses the block 
structure B • »B, approximates the eigenvector of the transfer matrix using a variational 
state in the tensor product 
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where £ 2 '''£iV-i denote m-state block spin variables. (Fig. la) The tensors A^* ^ and 

: ■< 

relation 



_B, , 3 are dependent on their positions % and j, and each of them satisfies the orthogonal 



where we have written s 1 and s N as ^ and respectively. Normally, they impose the 
normalization jl9| 

'V> t M+1 ) 2 = 1. (4) 
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The standard finite system algorithm improves the variational state (Eq.2) so that 
(Fig.lb) 

A (M) = J- vM, N TtfZVM N (5) 



all indices 



f 



is maximized via the tuning of Vt * , under the constraint Eq.4. The best Vt / 
is determined by the diagonalization of the renormalized transfer matrix 
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and by identifying its eigenvector Vt s s f with the tensor Vf M Sm+1 . The factor 

J <; Af-l s A/ s AJ + l>A/ + 2 «Af-l«Af + 2 

-l^A 



L e 1 and i? M+1 M+2 represent the renormalized half-row transfer matrix for the left and 



M + l 
M+l^M + 2 



the right half of the system, respectively. (Fig. lc) A pair of tensors A P M t and £>, * 

?A/-l^Af ?Af+l^ 

are then improved indirectly by rewriting the improved tensor M ^ M+1 using the singular 
value decomposition (SVD) 

2m 

yS M S M+1 = AS O B S M +1 (7) 

and by restricting the degree of freedom of £ M and £/vf+i down to m. The matrix £7^ ^ 
is a 2m-dimensional diagonal matrix p0| 



n = 



V 
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where the diagonal elements are in the decreasing order \u>i\ > |w 2 | > ■••|w 2 ml- The 
matrix Q satisfies the normalization Trfi 2 = J^f™ w | = 1- The finite system algorithm 
improves other tensors in Eq.2 by shifting the position of V by use of the wave function 
renormalization. pi], ^2| 

In the above standard improvement process for the variational state V s ■ M \ , the tensors 



A/ 1 t and Bt t are improved indirectly, only through the tuning for Vt t 

^M-l^M ^M + l^M + 2 ^A/-lV 

As a result, A P M p and B e P M+1 arc determined under the condition where 2m degrees 
of freedom is allowed for both £ M and Cm+i> although only m states are allowed for other 
block spin variables £ 2 • • ■ S, M _ 1 and £ M +2 " ' £jv- ^ ^ apparent that additional m numbers 
of freedom is allowed at the position where V is. (This excess freedom is common to both the 
finite and the infinite system algorithms.) Thus the variational state s is dependent 
on M, even after many sweeps of the finite system process. Figure 2 shows ln(A( M )) of the 
square lattice Ising model of the width N — 200 with free boundary condition, where we 
define the local Boltzmann weight as 

<fc (9) 
= «p{f(-i^ 1 + Vm+*i-J + ^ 1 -4)} 

for the Ising spins s = ±1. We keep (m =)8 states for the block spins. |^3| We have chosen 
the critical temperature K — J /k-g,T c . Since 2m degrees of freedom is allowed for both 

~ £/ s' s' 

£ M and £ Ar , q , the eigenvalue A^-* of the renormalized transfer matrix T ' M_1 M M+1 M + 2 [ s 
dependent on M; \W takes its maximum when M = N/2. 

The M dependence of the variational state V« a causes an ambiguity for the obser- 
vation of local quantities. For example, they calculate the local magnetization of the Ising 
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model using the formulation 



M = E v s [ m \ n s m vM n (io) 



and therefore (s M ) and (s M , ) for M ^ M' are calculated for different variational states 

V^ M g and Vs f 1 } N , respectively. The way to avoid such an ambiguity is simply to obtain a 
variational state that is independent on M. 

Now we show that we can further improve the variational state using the block structure 
B • B, and that the improved variational state is not dependent on M. The block structure 
B • B is known from the establishment of DMRG, jj], but the difference between B • »B 
and B • B has not been investigated from the view point of the numerical precision. In this 
case, the renormalized transfer matrix is constructed as 

rji^M — 1 S M*M+1 j^M -I s M ££ S m£m + 1 (11) 

where the improvement for the variational state is performed via the diagonalization of this 
2m 2 -dimcnsional matrix. The eigenvector U f , of the transfer matrix (Eq.ll) is 

2m 2 -dimensional, and it is possible to rewrite it as p2[ 

m 
m 

= yn' £ b £ " M (12) 

using the Gramm-Schmidt orthogonalization; Al M , and B f i M are orthogonal ma- 
trices that satisfy Eq.3, and Q'^ * and Q'^ , are m-dimensional diagonal matrices. 
(Note that in general fit t is not equal to Ql f .) Therefore, to use the block struc- 
turc B • B in the finite system algorithm is equivalent to write down the variational state 
using the tensor product (Fig. 3) 

at most m 

ui M \ = V A^, ...At" £ fit * (13) 
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and to improve A* M , and B f S , M directly via diagonalization of Tfw-i«UW < A11 

block spins £ 2 " " Cjv-i are a t most m-state; this is the non-negligible difference between 

U^\ N in Eq.13 and V S [ M \ N in Eq.2. 

Since £ 2 ' ' ' £jv-i are less than m-state in Usf* \ N , (Eq.13) we keep all the non-zero eigen- 
values of the density matrix during the numerical calculation of the finite system algorithm; 
we don't cut off any states. As we repeat the improvement for the tensors A^ £ and 

^ 3 for all i and j, the variational state S]v gradually approaches to its fixed point 

U s s , which is not dependent on M. As a result, unlike A^ A/ ' in Eq.5, the expectation 
value 

A' = V U s > a > T? 1 -?"!!, s (14) 

all indices 
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V u e s , e ff M - lS }j' M+1 u ( s £ 

^— ' -l" M^M + l Sm-I^M^M + I '■M-TM^M + l 

all indices 

is not dependent on M. Figure 4 shows In (A' ) of the Ising model; we use the same system 
size and parameter as those in Fig. 2. As we have expected, A' calculated from Eq.14 does 
not show any M dependence. It should be noted that A' is always larger than A^ M ) , which is 
shown for comparison. Since both A^ M ) and A' are variational lower bound for the partition 
function per row (= unit transfer), it is apparent that A' is better than A^ M -* in this case. 

We have used the block structure B • B to obtain A' after we calculated A< M > using 
the standard block structure B • »B. [^5| The additional calculation for B • B is not time 

consuming at all, because the matrix dimension of T e 1 * + in Eq.ll is smaller than 
rp i M-i m M+i m+2 ^ n Eq.6, and the finite system process can be done very rapidly with the 



use of the wave function renormalization. 21 , 22J Thus, treating B • B, one can increase the 



numerical precision in DMRG without increasing computational time so much. 

One might insist that the difference between Vs M \ in Eq.2 and U a „ in Eq.13 is 

1- ■ ■ N 1- .- JV 

negligible if m is sufficiently large. The statement is correct. We have to keep in mind, 
however, that occasionally people are suffering from increasing m in order to keep numerical 
precision. It is worth improving DMRG without increasing m, |^6| because to keep large m 
is difficult for complicated models. 

It is straightforward to introduce the B • B structure to the infinite system algorithm in 
order to recover the translational invariance of the variational state in the thermodynamic 
limit. The way is, as we have done for the finite system algorithm, first to employ the 
block structure B • »B and to perform the infinite RG processes till the renormalized wave 



function converges, |27 and then to use use B»B to perform additional infinite processes. In 



this way, the translationally invariant variational state by Ostlund-Rommer 15 is obtained 
numerically via DMRG. 
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tion, and I. Peschel for valuable discussion at the Max-Planck-Institute fur Physik komplexer 
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Figure Captions 

Fig.l Graphical representation of (a) the variational state, (Eq.2) (b) the variational eigen- 
value of the transfer matrix, (Eq.5) and (c) the renormalized transfer matrix. (Eq.6) 
Circles and squares denote bare and renormalized spin variables, respectively. We use 
black marks when the variables arc summed up in the corresponding equations. 

Fig. 2 The logarithm of \( M \ which is the eigenvalue of the renormalized transfer matrix 
ff M - lS ' MS ' M+1 f M + 2 of the Ising model at the critical temperature when N = 200 and 
m = 8. Open boundary condition is chosen. 

Fig. 3 The graphical representation of the variational state Us^\ N defined in Eq.13. The 

finite system algorithm for B • B gradually improves the state, and finally Us^\ N lose 
the M dependence. 

Fig. 4 The logarithm of A', which is the eigenvalue of T ' M ~ x M M+1 for the Ising model at 

?M-1 S M?M + 1 

the critical temperature when N = 200 and m = 8. Unlike ln(A^ M ^) shown for 
comparison, In (A' ) shown by circles does not show any M dependence. 
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